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Abstract 

The entanglement at finite temperatures are analyzed by using ther- 
^ ' mal models for colored quarks making up the hadron physical states. 

\Q , We have found that these quantum correlations entirely vanishes at 

T c > m 9 /ln(1.5). For temperatures larger than T c the correlations are 
' classical. Also we worked out the entanglement for the transverse Ising 

spin chain. In dependence on both temperature T and transverse field A 
we can identify a certain region, where the quantum effects are likely to 
dominate the system. We suggest the mutual information as a quantita- 
tive measure for the correlations in ground state. 
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p ■ 1 Introduction 

, One of the consequences of existing finite entropy at zero temperature is the 

understanding of Gibbs paradox, which represents one of very old thermody- 
namical problems. We now know that one has to take into consider the distin- 
' guishability for better understanding the finite entropy in ground state. On one 

hand side, we find that the change in the entropy, when non-identical particles 
have been mixed at fixed particle number and volume, goes to zero for vanish- 
ing temperature. This is known as the Nernst's heat theorem. On the other 
hand, for identical particles the entropy goes to a temperature-independent 
value 0(n) proportional to the number of mixed states at zero temperature. 
It is obviously independent on all thermodynamical quantities other than the 
number of particles itself. This is usually stated under the name of the third 
low of thermodynamics and dates back to the beginning of last century. For a 
non-comprehensible reason one used to mix up these two cases and suppose that 
in all systems the entropy goes to zero for T — > 0. This might be well founded 
from the accommodativeness to approximate S to zero for zero temperature. 

There are many physical systems, in which one has to avoid this approxima- 
tion. For instance, Schrodinger expected that the ground state of a gas of atoms 
encloses 2 N degenerate configurations in its structure, which must then provide 
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for it an entropy of AT In 2. The first application of these ideas on the quark 
matter has been introduced in pQ. It has been found that the SU(2) C color 
symmetry for each of the colored quarks in the ground state gives an entropy 
In 2. This has been extended to models at finite temperatures in Such an 
entropy can be seen as a reflection of the probabilities of quark mixing maneuver 
inside the mesons. Therefore, at zero temperature the confined quarks seem to 
be continuously tousled objects. 

A further example for the physical systems, in which a finite entropy exists 
at zero temperature is the recent lattice QCD results for the equation of state 
of quark- antiquark It has been found that the entropy takes the constant 
value 2 In 3 for T — ► 0. This is obviously well compatible with the evaluation 
of the entropy of SU(3) C colored singlet state based on the models given in |5]. 
There are other applications, for instance, the equation of state of the hadronic 
matter at very low temperatures 0], the quark-pair condensates at very high 
densities and low temperatures and the compact stellar objects [5]. 

In this paper we study the quantum correlations and their impacts in systems 
with SU (2) c color symmetry at finite temperatures. We suggest the concurrence 
as a measure for the thermal effects on these correlations. We aim for evaluating 
a temperature, up which the quantum correlations can survive. In section[3]we 
make a look at the transverse Ising model. The conclusions and outlook are 
given in section 



2 SU(2) C quarks states 

Using the models proposed in |2] based on the Boltzmann weighting for the 
finite temperature states of colored quarks we postulate that the total density 
matrix in the ground state of a meson singlet-state of SU(2) C color symmetry 
in the Hilbert-Schmidt space C 2 <g> C 2 could be written as 



PAB = j U®I-ae 't 5 (a z <g> I + 1 <g) a z ) - (l - e ^ ) ^ a t <g> a, J (1) 

I stands for the identity operator, e(p) = \J m 2 + p 2 is the single particle relative 
energy, a 1 ; i — x, y, z are the standard Pauli matrices and the quantity ae~ e ^ p ^ T 
describes the evolution of the color as a function of the temperature T. a 
represents the free parameter in this model. Clearly, for a — I we go back to the 
expression for the thermal reduced density matrix given in in equation (10). 
However, in order to fulfill the condition that the total density matrix has a 
non-negative value, we firstly suppose that 

"2^2 (2) 

For gl — 0, the mixed state in the whole system will be colorless for all 
temperatures. At T = equation ^ gives the expected density matrix for the 
pure colored singlet state in SU(2) C ground state. As introduced above, Eq. ^ 
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Fig. 1: Concurrence C as a function of T for two-site correlations in SU(2) 
quark states with m q — 5 MeV. For a — 0, the 'critical' temperature takes the 
value T ~ 12.5 MeV (dotted line). For a -> 1/2 we find that T c -> oo (solid 
line), meanwhile for a = 0.499 we get T c = T c d (dashed line). 

is also able to reproduce the finite temperature behavior, where we expect no 
quantum correlations to be exist. There is no direct description of a kind of 
quantum transition available by means of Eq. ^ One should notice that pab 
satisfy the condition Tr p 2 < 1. 

In order to compute the entanglement of such states at finite temperatures, 
we use the concurrence formalism which is defined as 

C = max{Ai - A 2 - A 3 - A 4 , 0} (3) 

where A's are the square roots of the eigenvalues in decreasing order of 
PAB{py ® a yPAB a y ® a y)- Pab i- s the corresponding complex conjugation in 

the computational basis {| + +), | H — ), | — h), | )}• After straightforward 

little algebra we find that 

C 1-(1+ Vl^_4ai) e -^ 1 if T<T C 

C(T) = I - (4) 

[ if T > T c 

where eSm - Tc - 
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Fig. 2: t = C 2 as a function of the linear entropy Sl — 4/3(1 — Trp 2 ) for the 
Maximally Entangled Mixed States (MEMS) and for the entropy calculated for 
a = and 1/2. 



The value of a can then be fixed through the hypothetical assumption that 
the transition temperature from confined hadron gas to deconfined quark 
gluon phase at zero chemical potential would have the same value as that for 
the transition from entangled to non-entangled SU(2) C colored states. In this 
case 

« 2 = j(l-4( e ^-l) 2 ) (5) 

From the lattice QCD simulations we know that for two quark flavors the tran- 
sition temperature at /j, q = is ~ 173 MeV. If we take the mass threshold rn q (/j,) 
of light quarks as e(p) [B] where \x in this case is a renormalization group scale 
and for simplicity we set m q (p(T)) = 5 MeV, we therefore find that a for which 
T c = T c d is 0.499. 

As shortly given above, the entanglement transition is corresponding to C = 
0. Therefrom, we can simply estimate T c for the simple case that a = and 
m q = 5 MeV, as 

T c - 12.5 MeV < T c d (6) 

This represents the lower value of T c . In Fig. we draw C for the two-site 
correlations in SU(2) C quark states as a function of T with different values of 
a. We find that T c <C T r d for a = 0, meanwhile T c — > oo for a — > 0.5. 
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Fig. 3: Single-site color entropy as a function of the temperature T (see text 
for more details). 



The density matrices given in Eq. ^ have some particular properties. For 
specific values of T and a (see Fig. |2Jl they belong to the states with high 
entanglement, the so-called Maximally Entangled Mixed States (MEMS) UJ. 
This would suggest that the thermal evolution of the colored quark will maximize 
the entanglement at given degree of mixing. 

Moreover, it has been shown in that in an antiferromagnetic ring for 
Heisenberg model (one-dimensional with boundary conditions) with an even 
number of spin- 1/2 the maximum entanglement is equal to the entanglement of 
the ground state (minimum energy) with zero z component of the spin. The 
transverse spin models will be studied in section [3] Therefore, we do not need 
for the moment to elaborate any further details. 

In Fig. [3 we plot r = C 2 versus the linear entropy Sl = f (1 — Trp 2 ). 
Clearly, we notice that the states with the properties given in Eq. ^ typically 
do not maximize the entanglement at a given degree of state mixing. These are 
the states with a = 0. But for the special case a — > 0.5 such states are very 
close to maximize the entanglement at all degrees of mixing. For Sl — 0, i.e., 
pure state, the results for MEMS and all values of a are coincident. 

As shortly given in section^ the effects of the quantum entropy of SU(2) C 
of colored quark states on the bag pressure have been discussed in [3]. This 
we can also be analyzed here by computing the color entropy of single-quark 
in one meson. In Fig. we show the results for different values of a. In the 
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Fig. 4: The relative entropy of entanglement is depicted in dependence on 
T (a = and m q — 5 MeV) . The quantum correlations are given in Eq. 
The classical correlations are the difference between the total correlations from 
the von Neumann mutual information and the quantum correlations. At T c the 
quantum correlations are entirely vanishing meanwhile the classical ones remain 
finite for all temperatures. 

case a = we notice that the quark quantum entropy S q ,2 is a constant for 
all temperatures. As given in 0] this case is corresponding to excluding the 
quantum entropy from the equation of state. The reason for this result is clear, 
since for a = all states are composed of completely random swfetates mixed 
together with singlet states. In both groups the single-site entropy is maximum 
(the net color is zero). Though when a — there is no color polarization 
for each quark (local polarization is zero) and the overall color is also zero, 
the colors of single-quark are still correlated. For temperatures below T c the 
correlations is a mixture of quantum and classical ones. However, for T > 
T c the correlations are completely classical. To illustrate these features we 
plot in Fig. 0] the entanglement measured by the relative entropy from the 
entanglement, which is defined as 

E r = minS(p||p*), (7) 
For completeness the discussion, we mention that the classical correlations are 
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defined as the difference between the total correlations measured by the von Neu- 
mann mutual information and the quantum correlations measured by the rela- 
tive entropy from the entanglement |12| . 

*(p) = S(p\\p A ® pb) - min S(p\\p*), (8) 

where T> is a set of all separable states in the Hilbert space, in which the 
su&spaces A and B are also defined, pa and ps are the corresponding reduced 
density matrices, respectively. The upper limit of temperatures, at which both 
quantum and classical correlations exist, is also drawn in Fig. ^ as dotted line 
(a = 0). As expected, at T c the quantum correlations are zero meanwhile the 
classical ones remain finite. That is the case also for higher temperatures. 

3 Transverse Ising model 

In this section we would like to investigate the quantum correlations (entan- 
glement) for some spin models in which strong correlations exist and compare 
their results with the results of the SU(2) C colored quark states given in pre- 
vious section. The general category of all spin models goes under the name of 
the Heisenberg model. For transverse Ising model with N sites the Hamiltonian 
reads 

JV-l 

H = -J2(^v? + i+<), (9) 

i=0 

where a 1 , i = {x,y,z} are the standard Pauli matrices and A is the transverse 
field. It should not be confused with the eigenvalues given in Eq. |3J At 
zero temperature it is expected that this system goes through a quantum 
phase transition QPT (from para- to ferro-magnetic phase transition). A < 1 
corresponds to a para-magnetic and A > 1 to a ferromagnetic phase. Unlike the 
ordinary phase transition, which occurs at finite temperatures, the fluctuations 
in QPT are fully quantum and the transition is due to the changes in the 
ground state wave functions. At the critical point long-range correlations are 
highly expected. This kind of correlations is also expected near the transition 
point in the classical thermal phase transition. However, as the system has zero 
temperature and with the assumption that the ground state is non-degenerate, 
the system, in whole, is expected to build up a pure state. It follows that 
the correlations, which arc in principal the accessible signatures for QPT, are 
obviously stemming from the long-range entanglement in the ground state. 
Through the change in the correlation functions which will be reflected in 
fundamental changes in the entanglement in ground state, we get a distinct 
signature for probing QPT. 

Before we study the two-site entanglement we want to know the role of 
A in describing the entanglement. To illustrate this we start in Fig. [S] with 
the single-spin quantum entropy which in turn brings about a device for 
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Fig. 5: Zero temperature single-site entanglement S of the ground state of 
transverse Ising model (solid line). The dashed line gives the corresponding 
results for the thermal ground state. 

measuring the entanglement of one spin with the rest of the chain. This can 
be given by evaluating the entropy of the single-site density matrix pa of the 
ground state after tracing out all other spins (see appendix) . In transverse Ising 
model we notice that S varies from zero at A = to a maximum value at A = 1 
(Fig.[SJ). For higher values of A, S declines, rapidly, and for A — > oo, S reaches 
the asymptotic value zero. For the thermal state (pi = limr_»o Tr^e^^/Z) 
we can also apply the von Neumann prescription for calculating the entropy. 
We get here the same behavior as in the ground state entropy for < A < 1. 
But for A > 1, we notice that S goes on with the increasing until it reaches the 
asymptotic value 1 for A — * oo. This is because of the realization of equal mixing 
of spin thermal states. We can so far conclude that the single-spin entropy in the 
thermal ground state is not measuring the entanglement in the limit A — > oo. It 
merely gives the degree of mixing. For the first time, we study here the single- 
site thermal quantum entropy and its derivative with respect to A. It should 
be understood that in this case the quantum entropy is not a measure for the 
entanglement, however one can expect that the change in the quantum entropy 
could, in turn, reflect the change in the entanglement. In Figs. E] and we 
show the results in kT — A plan. We notice that when A = 1 both functions are 
approximately maximum (VT). As expected, the first derivative of the quantum 
entropy for high temperatures is always picked at A = 1. This region is not fully 
compatible with the thermal concurrence published in ,13 . Concretely, when 
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Fig. 6: Single-site thermal entanglement S as a function of T and A in transverse 
Ising model. 




Fig. 7: First derivative of single-site thermal entanglement d\S as a function 
of T and A. 




Fig. 8: The entanglement measured by the concurrence of nearest- and next- 
nearest-neighbor two spins at zero temperature depicted in dependence on the 
transverse field A. The next-nearest-neighbor concurrence is very small. Here 
it is multiplied by factor 10. 



we compare our results (Figs. |H] and 0) with the Fig. 6 in I n the latter we 
see that for high temperatures the concurrence increases for decreasing A! 

The two-site entanglement, as its name indicates, measures how two spins 
separated by a distance r are entangled. This measure is to be accomplished by 
evaluating the concurrence of the entanglement of formation of two-site density 
matrix after tracing out all other spins in the chain (see appendix). The two- 
site entanglement for nearest-neighbor and next-nearest-neighbor spin-pairs at 
zero temperature for the transverse Ising model are given in Fig. |H1 The length 
of entanglement £ aparently vanishes for all other pairs, since the concurrence 
gets negative. Therefore, the true non-local quantum part is a short-range one, 
meanwhile the classical part is of long length (diverging for an infinite system) . 
We now would like to compare the classical with the quantum correlations. 
We use the same formalism as in the previous section, however, since it is 
not possible to directly evaluate the relative entropy for the two-site density 
matrix, instead of this we can bound it (upper with lower bounds) to estimate 
the quantum correlations and afterwards extract the classical ones. The upper 
bounds are calculated from the entanglement of formation 

j^^( (1 + V p^-) do) 
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Fig. 9: Bounds on classical and quantum correlations for transverse Ising model 
are drawn versus the strength of the transverse magnetic field A. 



where h(x) — —xlogx — (1 — x) log(l — x) is the Shannon entropy. The lower 
bounds read 

E?» = (C - 2) In (l - |) + (1 - C) log(l - C) (11) 

In Fig. the bounds on classical and quantum correlations are depicted in 
dependence on A. We see that the quantum correlations are smaller than the 
classical ones. The latter are merely produced by the entanglement in the 
ground state. This suggests the study of the mutual information as well as 
the single- and two-site entanglement [161 117) . because the mutual information 
combines all correlations (classical and quantum) which are created by the 
entanglement |17| . 

In ^H] the scaling of concurrence has been studied. It has been found that 
for infinite system the first derivative of the concurrence d\C is a diverging 
function at the critical point. It is important to see whether the mutual 
information has such a property. We show in Fig. the first derivative of the 
mutual information d\I, top curve. We obviously get for the first derivative 
of concurrence d\C, bottom curve, the same dependence as we obtained for 
d\I. This has been studied elsewhere ^Hj- As expected, at the critical point 
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Fig. 10: First derivative of the quantum and classical correlations as functions 
of A. For A < 1 both derivatives are finite. But they are vanishing for A — > oo. 
At the critical point (A = 1) both functions are divergent, indicating radical 
change in the ground state structure. 

A = 1 we notice that d\I is singular and divergent. This divergence is a 
non-ambiguous signal for the change in the ground state structure. Therefore, 
we expect that the mutual information can also be used, in order to probe the 
entanglement in ground state. 

Thought the color density matrix Eq. ^ as discussed in the previous section, 
does not undergo QPT, it is still possible to be compared with the transverse 
Ising model at zero temperature. In fact for A > 1 the mixing due to the change 
of the temperature in colored system is compatible with the degeneracy mixing 
in the ground state. Clearly, for A > 1 the relative entropy (Fig. 21 IH1 an d EJl 
in both systems are decreasing functions for increasing T and A. However, the 
classical correlations are different in both systems (see Fig. 01 and EJ. This is 
due to the tendency of the color system to be in a completely random state by 
loosing all kind of correlations at high temperatures. 

4 Conclusions and outlook 

In the limit of the proposed thermal models for colored quarks making up the 
hadron physical states we have utilized the concurrence, in order to study the 
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quantum correlations in SU(2) C colored quark states. The aim of this study is 
to see whether the quantum and/or the classical correlations can survive the 
deconfinenient phase transition from hadron to quark-gluon plasma. Obviously, 
we have found a region, where an entanglement transition is highly expected. 
This gives an upper temperature limit T c , below which the quantum correla- 
tions are significant. At temperatures higher than T c the quantum correlations 
entirely disappear allowing dominating classical correlations. 

We have also studied the correlations in another system, that has the same 
SU (2) symmetry, namely the transverse Ising model. Using the single-site en- 
tropy we have identified a certain region in the T — A plan, where the quantum 
effects likely dominate the behavior of this spin system. The classical correla- 
tions can be extracted from the knowledge of quantum and total correlations. 
We can conclude that the nature of quantum correlations is short-ranged, mean- 
while that of classical correlations is long-ranged. We have also shown that the 
mutual information, which combines both types of correlations, is a possible can- 
didate to study the structure of QPT. In a forthcoming work we shall address 
the mutual information and its impact in describing the corrlations at finite 
temperatures in more details. Also we shall utilize the other models given in 
for the SU(3) C quark states, which aparently is more complicated than SU(2) C . 
It might be possible to compare the quantum phase trnasition in SU (3) c colored 
quark states with the possible transition in the interior of compact stellar. 
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5 Appendix 



For the one-dimensional quantum Ising model with N spin-1/2 sites the one- 
and two-point correlation functions cab be calculated by using the operator 
expansion for the total density matrices in terms of the Pauli matrices. 



Pk,i = 



fc=0 

1 3 

7 I] (^W>°-fc®°; 



(12) 
(13) 



i=0,j=0 



As given in |14| we can replace the sums by integrals and afterwards calculate 
the different correlations in dependence on the separations between arbitrary 
two-sites, r = \k — l\. Then for the thermodynamic limit N — > oo we find that 



(of) = -- I dk(l + Acos(fc)) 
where /3 = 1/T. For T = we have 

«) = 

Using the definition: 



tanh(/3w/2) 



0, A < 1 

( 1 _ A -2)l/8 j A>1 



(14) 



(15) 



- / dfccos(fcr)(l + Acos(fc))- 



— — I dk sin(fcr) sin(fc) 



tanh(/3w/2) 



(16) 



then the two-site correlation functions read ^1 fT^j 



T yrrV\ - 



C-l C-2 

Co C-i 

C r -2 C r _3 

C2 C\ 
Cf Cf — 1 

7 Cl) 2 ~ C r C— r 



C — f 

C-r+1 

. C-x 

C-r+2 
C-r+3 



(17) 



(18) 



(19) 



It should be noticed that in Eq. ^5] there is no factor 4 in the front of the first 
term as the case in As in ^3] we find that the factor 4 should be multiplied 
by the magnetization m z and not by the polarization (<x z ). Since m z — (a z )/2, 
we get factor 1. 
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